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I. INTRODUCTION 



The soft-collinear effective theory (SCET) has been proposed systematic 

framework for the study of processes involving energetic hght quarks and gluons. Possible 
applications include the decays of heavy hadrons into light particles in the kinematical 
regions where the final products are very energetic, and hard scattering processes involving 
light hadrons, such as deep inelastic scattering and exclusive hadron form factors at large 
momentum transfer. SCET provides a natural framework for establishing a systematic 
expansion in A/Q where A is the QCD scale and Q is the typical large energy of the particles 
involved. In particular, it provides a convenient tool to establish factorization theorems and 
study power corrections. 

The observation that lies at the basis of SCET is that all the kinematical singularities 
appearing in these processes are connected to the exchange of collinear and soft particles 
("long distance modes"). Therefore, an explicit description of these processes in terms of 
these degrees of freedom offers the usual advantages of an effective theory approach. On the 
one hand, the power counting is greatly simplified since it can be performed at the operator 
level; on the other hand, standard renormalization group techniques can be easily applied. 

For definiteness we consider the decay of a heavy quark of mass m which emits light 
quarks and gluons of energy Q c:^ m. The possible external states are the soft spectator 
(whose momentum is of order A) and the decay products (light mesons and leptons), in 
addition to the heavy quark which cannot give rise to singularities and can be considered 
as an external source. By the Coleman- Norton theorem j5|, the infra-red singularities of 
Feynman diagrams describing these amplitudes correspond to the propagation of on-shell 
particles. The idea underlying SCET is to identify all the possible on-shell modes that can 
appear in the initial and final state (and that are, hence, responsible for all the infra-red 
singularities) and to write the effective theory of their interactions as an expansion in A/Q. 

In the following we use the standard light-cone decomposition of momenta 

1 _ 1 _ 

= -n^n ■ p + -n'^n ■ p + p'^ = {p+,P-,p±) , (1) 

where n and n are light-cone vectors satisfying n'^ = n"^ = 0, n ■ n = 2. In any given process 
n and n are chosen to be aligned to the final state collinear momenta. We also introduce the 
dimensionless parameter A that will serve as the expansion parameter of SCET. Its precise 
definition is discussed below. 

In the description of heavy meson decays, there are three relevant kinematical configura- 
tions. 

a) Soft quarks {qs,hy*) and gluons (A^) with momenta Ps — A = Q{X,X,X) (where 
we defined A = A/Q). The propagation of these modes is described by the non-perturbative 
regime of QCD {pi ~ A^) and is therefore incalculable. Therefore, the exchange of soft 
particles can only be parameterized and results in the non-factorizable contributions to 
heavy-to-heavy and heavy-to-light form factors and the B meson wave function. 

b) Collinear quarks (^„) and gluons (A^) with momenta Pc — (5(1,A^,A). These 
modes appear in the description of the constituents of a fast-moving light meson. They 



* hy is the usual heavy quark degree of freedom that, after removing the fast osciUating components, is 
equivalent to a soft field. 
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are also non-perturbative and their exchanges can only be parameterized in terms of the 
light-cone wave functions of the final state mesons. 

c) Hard collinear quarks (S„) and gluons (A'^^) with momenta p/^c — Q(1,A,V^)- 
These modes are necessary to describe inclusive jets {e.g. the Xg system near the end-point 
region of the photon spectrum in the inclusive B Xs'j decay jlf) and interactions of soft 
fields with collinear particles (e^. the 5- meson soft spectator after being struck by the 
energetic photon in B —>■ jeu ja 0, B 1^)- Hard collinear modes have virtuality of order 

~ AQ ^ A^, can be integrated out perturbatively and result in so-called jet-functions. 

Further intermediate 'hard' collinear modes with momenta Q{1, A, A) occur at tree level 
(for instance in the coupling of a soft and a collinear field). It is not necessary to introduce 
explicitely such fields, because they are connected to hard collinears via a reparameterization 
transformation. In Refs. ^ three different invariances under changes in the light cone 



vectors n and n were introduced: 

type I 
type II 
type III 



n 



" ^n'' + A^, ^ (2) 
nf^ ^nf^^ e^, ^ n^" (3) 
fif'/a (4) 



where ~ 0{X), (e^,a) ~ 0(1) and A refers to the expansion parameter relevant to the 
theory considered. In our case with two types of collinear fields, there are actually two ex- 
pansion parameters, ^JK/Q and K/Q for the collinear and hard collinear fields, respectively. 
Thus we must have two invariances, related to the corresponding fields. In particular, the 
hard collinear sector must satisfy a type-I reparameterization invariance with A^ ~ 0(-\/A). 
This invariance connects the (1, A, V^) and (1, A, A) modes. 

From a technical point of view, it is always possible to completely integrate out the hard 
collinear modes because in all applications they always appear as internal modes. This 
is obvious in processes like B —>■ '-feu. In more complicated situations like heavy-light 
semileptonic decays B — >■ neu, a two step procedure has been proposed In the first step 
an effective theory is formulated (SCET-I) containing only usoft and hard collinear modes, 
which is matched into a second step onto the final effective theory (SCET-II) containing the 
collinear and soft modes. [In the alternative treatment of Ref. [l3| only the modes a) and 
b) are introduced and one matches directly from QCD onto SCET-IL] On the other hand, 
in inclusive decays (like B Xsj) one can always write the decay width using the optical 
theorem and again hard collinear modes can be viewed as internal particles. 

For the purpose of this paper, we take the point of view that problems associated with 
the integration over the hard collinear modes at subleading order have been cleared and 
proceed to the analysis of power suppressed contributions. Once the process is specified, 
the integration of the hard collinears gives, order by order in a<j, all the relevant SCET 
operators (that will involve only soft and collinear modes). The matrix elements of these 
operators between initial and final states fall in two groups: those factorizable in terms 
of "conventional" form factors and light-cone wave functions, and others that require the 
introduction of new non-perturbative objects. In particular, the light-cone wave functions 
enter through matrix elements of SCET operators (involving two quarks and an infinite 
number of gluons) between the vacuum and a meson state. 

In order to trust this perturbative computation, one also has to show that higher loop 
contributions will not introduce new non-perturbative structures at a given order in A. 
This can be achieved by writing the most general set of operators allowed by gauge and 
reparameterization invariance and showing that their matrix elements factorise. 
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In this paper we classify all the possible SCET coUinear operators that appear at leading 
and subleading order in A and compute their matrix elements in terms of the usual light 
cone distribution amphtudes of pseudoscalar ji^, 21 1 and vector 2^ 2^ 2^ mesons. These 
operators are necessary for the SCET analysis of any process involving energetic light mesons. 
In particular, we find that for decays involving transverse polarized vector mesons, only 
subleading operators contribute and that, in this case, it is not possible to neglect the 
contribution of three-particle distribution amplitudes. 

Note, finally, that there are two equivalent formulations of SCET, that are usually denoted 
as hybrid [l|, and coordinate space [13, Il3 5 respectively. In the former, strongly oscillating 
collinear modes are removed from the theory by a partial Fourier transformation; in the latter 
the slowly varying soft fields are multipole expanded. These manipulations are necessary in 
order to fully expand the Lagrangian in powers of A. A failure in achieving correctly the 
complete A-expansion would result in a theory whose infrared behaviour does not reproduce 
the one of full QCD j2^. We will work in the hybrid formalism and show how to readily 
translate our results into coordinate-space. 

We start by summarising briefly the main ingredients of the hybrid formulation of SCET. 
In Sec. ini we present the complete list of leading and subleading operators built only of 
collinear fields. In Sec. IIIII and IIVI we show how to extract the matrix elements of the 
various SCET operators and present our results. In Sec. |3 we show that reparameteriza- 
tion invariance (RPI) gives strong constraints on the Wilson coefficients of the subleading 
collinear operators. In Sec. IVII we present a sample application for the subleading collinear 
operators constructed in this paper, discussing final state photon emission in weak annihila- 
tion. We check our findings for the matrix elements of these operators by verifying an exact 
Ward identity; this illustrates the numerical impact of the three-particle contributions which 
are usually neglected in practical computations. We collect our conventions in Appendix 1X1 
and in Appendix ^ we give the parameters of a minimal set of light cone wave functions 
satisfying the equations of motion in QCD used in the numerical evaluations in Sec. IVII 



II. SUBLEADING COLLINEAR OPERATORS IN SCET 

In this Section, we present a complete set of collinear quark operators ^ri ■ ■ ■ ^„ up to sub- 
leading order in A. Such operators are important in their own right, as interpolating fields 
for meson states, and as building blocks for more complicated operators in the effective 
theory. In constructing a complete basis of operators we use constraints from the collinear 
gauge invariance, Dirac structure of the collinear quark fields ^n, and from reparameteriza- 
tion invariance of the effective theory. A similar procedure was recently used in Ref. to 
construct the complete set of heavy-to-light currents in SCET. 

The operators to be constructed in this Section contain the collinear quark ^n,p and gluon 
An^q fields, together with the collinear covariant derivative iD^ = + gA':^. In order to be 
able to perform the power counting of the various contributions at the operator level, it is 
convenient to assign a A scaling to the fields requiring the kinetic terms to be of order 0(1). 
In this way, one obtains ~ A and A'^ ~ (A^, 1, A). 

In the hybrid formulation of SCET, one achieves a correct A expansion by extracting the 
large Fourier modes from each field: 

0,(x) = 5^e-^^^>,,,-,(x) (5) 
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where pc = Q{0, 1, A), are labels and the new field 0c,pc is responsible for fiuctuations with 
momenta of order Q{X^, A^, A^). 

It is convenient to introduce a "label" operator |3|] which acts on the collinear fields 
and picks up their large momentum: ^n,p = {■^n^ + P±) ^n,p and V^, respectively. When 
acting on a product of several fields, these operators give the difference between the total 
label carried by the fields minus the total label of the complex conjugated fields. 

We will also use a special notation which associates a momentum label index to an 
arbitrary product of collinear fields. Our convention is 



Xn,i 



[W^^^l = [S{uJ - n-V)W^^ , [WHD^^W]^ = [6{uj - n-V)WHD^,W] (6) 



where 6{uj — n ■ V) acts only inside the square brackets. 

The collinear operators that we write in the following can include a nontrivial fiavor 
structure. When required, this will be denoted by a superscript showing the quark fiavours. 
Note, finally, that for each operator we can have the singlet and octet colour structure. We 
will write explicitly only the former since the matrix elements of any octet operator between 
the vacuum and a meson state vanishes. 

At leading order in A there are only three independent collinear operators, which can be 
chosen as 



T ( -'\ — I 



i. 



(7) 
(8) 
(9) 



where 7" = 7" — n"^/2 — n"^/2 and uj = [001,002)- Their transformation properties under 
charge conjugation are 



J'yj){uJl, OO2) 



'^VT (~"^2j —^l) 



(du). 



(10) 

(11) 



At subleading order in A, the number of allowed structures is much larger. It is convenient 
to choose a basis of collinear operators with simple transformation properties under charge 
conjugation. We choose the following four chiral-even collinear operators 





'-4 




V^O) = 


'-4 




AU^) = 


4 






'4 





1 



n-V 



X 



n-V 



U)2 



1 _ 1 

7°75Xn,<^2 + X n,u>-i 



1 _ 1 



n-V 



n 



n-V 



(12) 
(13) 
(14) 
(15) 
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and the three chiral-odd operators 



n 



n-V 



75 ^ 



/9 /3 
7_L _ - _7_L 

rXn,(^2 Xn,uji _ 



.-pt- 



ljJ2 



(16) 
(17) 
(18) 



U2 



together with the corresponding colour octet operators which will be denoted by the same 
letter and a colour index a. We include the factors 1/n-V and in the definition of 

the operators, Eqs. ()12|) - ()18|) . to make them invariant under the transformation n na, 
n ^ n/a (type-III reparameterization invariance). 

These operators are not the most general coUinear gauge invariants at 0(A). In analogy 
to the heavy-to-light current considered in Ref. [ll|, it is possible to write also three-particle 
operators, which contain three coUinear gauge invariant factors. Their Dirac structure is 
again restricted by the effective theory constraint iji^n = 0, which leaves two possible chiral- 
even operators 



n,u}2 Xn,uii 



n-V 



WkgnaGf^'^W 



X' 



n,uJ2 



(19) 



Xn,Lui n 75 [ _ ^ 

2 n-V 



1 "A 



n-V 



and a single chiral-odd operator 

rf{0) = I 7! [^W^tDiWU Xn,.2 = Xn,.r\ll 



n-V 



Xn,ui2 

(20) 

Xn,uJ2 - 

(21) 



The factor 1/n-V assures again invariance under type-III reparameterization invariance. 
The transformation properties of the subleading operators under charge conjugation are 



5("^)"(^i,^2) 

P("'^)"(^l,^2) 



-VI2' {-UJ2, -UJi) , 

^1,2 {-UJ2,-UJi), 

5W"(_^2,-^l), 
_pW«(_^2,-^l), 

TW"/5(-cu2,-u;i). 



(22) 
(23) 
(24) 
(25) 
(26) 



The corresponding operators of opposite charge conjugation properties can be constructed 
by changing the relative sign of the two terms in Eqs. (fT ^ -(fTH j) . They will be denoted with 
a tilde, e.g. 



vr(^) 



7 Xn,uJ2 Xn,uji1 _ 



n-V 



(27) 
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and similarly for the remaining 6 operators. They are not independent and can be related to 
the above ones by using the Dirac identities in the Appendix A. For example, using Eq. ()A6|) 
(written in terms of the transverse antisymmetric tensor = ^e^uaisn'^n'^) 

f7M = ^?r|-<^|75, (28) 

one obtains the relations 

Vf = + le^^Mu , At = A'i + le^lV^u , (29) 
= Vl^ + ie^lA2u , = + zef V2. . 

They can be solved for Vi, Ai, with the results 

VJ=iefAip, AJ = iefVip, (30) 

= ief{Aip - A2p) , = lefiVip - V^p) • (31) 

For the chiral-odd operators there are two identities following from their definition 

9i,r'^'{^) = s{^) , g^.f'^'i^) = s{^) , (32) 

and two other identities following from Eq. ()A8jl 

V = 2ei,r^^ , V = zei^T^-' . (33) 

Let us finally discuss the way these operators appear in explicit calculations. A given 
QCD operator Oqcd is matched onto SCET operators containing the subleading collinear 
bilinears introduced above 



Cqcd — ■ ■ ■ + 



j duJiduj2Ci{ui,U2){- ■ ■}Vi{uJi,uj2) (34) 



+ J duJidLJ2duJ3C2{uJl,UJ2,UJ3){- ■ ■}Vi{uJi,U2,UJ3) 

where the ellipses {■ ■ ■ } denote possible soft fields which were omitted in writing the SCET 
operator. The Wilson coefficients Ci_2(cUj) depend on the momentum labels of the collinear 
bilinears. 

After factorization, the matrix elements of the collinear operators Vj(co'j) between a light 
meson and vacuum lead to non-perturbative functions {M{pM)\Vi{ui,uj2)\0) ~ (pi{u). It is 
convenient to implement momentum conservation uJi — lj2 = Pm by introducing the mo- 
mentum fraction u by (ti;i,ci;2) = {u,—u)n ■ pm, with u = (0,1). The charge-conjugation 
transformation properties of the collinear operators Vi(ti;i, ti;2) (see Eqs. flTm) . (fTT|) . (P^ - (Pm) ). 
taken together with the C quantum number of the state \M{pm)), fixes the symmetry prop- 
erty of matrix elements under the substitution u ^ u. For example, taking C = — 1 as 
appropriate for the p meson, one has 

(P(P, r^Wr"-"''^ (^1, 0^2) |0) ~ ^f"^--^ (u) , (35) 

such that only the odd (even) part of the corresponding Wilson coefficient C{lji,U2) will 
give a non- vanishing contribution to the given matrix element of Eq. (jMj) . The matrix 
elements of Eq. fl34|) will be given with the integration measure Iljda;, replaced with du for 
the 2-parton operators, and with dQ!idQ!2dQ!35(l — E^^^^aj) for the 3-parton operators. 
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III. STRATEGY FOR THE COMPUTATION OF THE MATRIX ELEMENTS 



In this section we explain the procedure that we use to extract the matrix elements 
between the vacuum and a state with one vector or pseudo-scalar meson of the leading and 
subleading collinear operators introduced in Sec. [TTl We will see that they can be expressed 
in terms of the vector and pseudo-scalar mesons' light-cone wave functions. 

The usual starting point are matrix elements of nonlocal operators of the form 

{M\q{x)VW[xMq{y)\^) or {M\q{x)VW[x,z]G^,{z)W[z,y]q{y)\Q). (36) 
The W[-, ■] terms are Wilson lines required by gauge invariance and are explicitly given by 



W[x, y] = Pexp 



ig{x-yY / dtA^{tx + ty) 



(37) 



where t = 1 — t. In the following, the presence of the appropriate Wilson lines between fields 
evaluated at different space-time points is understood. 

The general procedure that we adopt is rather simple. We first project the QCD oper- 
ators onto SCET keeping leading and subleading contributions; then we use the standard 
definitions of the light-cone wave functions given in Refs. HQHIilliJ to identify the 
matrix elements of the various SCET operators. 

The presence of Wilson lines and of fields situated at different space-time positions, causes 
the appearance of arbitrarily suppressed operators in the projection of the QCD operators 
onto SCET. In order to simplify the extraction of the leading and subleading operators we 
expand the QCD operators around the transverse direction: 



q{x)Tq{y) 



q{x)Tq{y) 
d 



d 



X I 



q{x)Tq{y) 



dx'' 



yj_=o 



q{x)Tq{y) 

y±=o 

f higher powers of x±, y± , 



where 



d 
dx'^ 
d 



q{x)Tq{y) = q{x) D ^^T q{y) - i / dit z"" q{x)gGo,^{tx + ty)T q{y) , 



dy[ 



jrq{x)Tq{y) 



q{x)T D±^q{y) - i dt t z"' q{x)gGa,^ 
Jo 



{tx + ty)Tq{y) 



(38) 

(39) 
(40) 



and we define z = x — y. 

All the terms of this expansion count as 0(1) (in fact, we only require x — ?/ to be on the 
light-cone without any constraint on the direction of this vector) but the projections of the 
various terms start with SCET operators of increasing A suppression: 



d 



dx^^ 
d 



dx, dx\ 



q{x) 'W[x,y]Tq{y) 
q{x)W[x,y]V q{y) 
q{x)W[x,y]V q{y) 



!/±=0 



00 + 01 + ■ ■ 

01 + O2 + ■ ■ 

02 + 0^ + ■■ 



(41) 
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where On are generic SCET operators suppressed by a factor A" compared to Oq. Since we 
are interested in leading and subleading operators only, we just need to consider the few 
terms explicitly written in Eq. (jHHj) . 

Let us now show in detail how to perform the SCET projection of a given QCD operator. 
The starting point is to express the QCD quark field in terms of coUinear and soft fields 



Ps 



Pc 



in ■ D, 

Since all the terms in the expansion ()38p are evaluated at x_l = we obtain 



where we do not include the soft quark because, as argued below, its contribution receives an 
additional power suppression and it is negligible for the analysis of leading and subleading 
operators. 

Let us now consider the leading order projection of the QCD operator q{x) W[x, y\ 7^ q{y)- 
The matrix element of the term is given by 

(M(pm)I ^V^uJXn,^,n^\xn..M = ^'d«e^("P«-^+^P--^)(M|x„,„p,,n'^|xn,-.,„|0)(44) 

oc /" dMe'("P"-"+"P"-^)(^(M), 
Jo 

where 

''v^LU= /"dcuidcj2e^("i"™"-^), (45) 



and u = 1 — u. The Wn factors (contained in the Xn fields) come from the projection of 
the full QCD Wilson line, Eq. ljTfjl . in the limit x± = = y±. In Eq. (jllj), we replaced 
the discrete sum of Eq. ()43|1 by the integral over Ui. Once momentum conservation is 
imposed, the sum of the two labels must be equal to the total momentum of the meson: 
uji — UJ2 = n ■ Pm and the fractions of momentum carried by the two collinear quarks are 
^1 = upM and U2 = ~upM, respectively (this assumes that the parton carrying momentum 
LJi{iJ2) ends up as a quark (antiquark) in M). The equality in Eq. (jl^ is obtained by direct 
comparison with the definition of the QCD light-cone wave functions for the q'j'^q current. 
The result is 

(M|x„,^,n''|xn,..jO) (X <f{u). (46) 

Note that although in principle mixed soft-coUinear terms can be generated too, their 
matrix elements on a collinear meson state require insertions of the soft-collinear subleading 
Lagrangian and are suppressed relative to those of the diagonal operators we consider 
here. 

Before concluding this section, let us explain how to translate the above arguments into 
the coordinate-space formulation of SCET. In the first place note that our leading and sub- 
leading operators do not involve soft covariant derivatives. This implies that each operator 



9 



can be trivially translated into coordinate-space formalism by replacing the label operators 
with ordinary derivatives: V^+g A^^ id^+g A^^ and n-V in-d. The latter substitution 
might rise concerns about the invariance of the coordinate-space version of Eqs. (fT^ -()21 |) 
under soft gauge transformations l2| • From the discussion in Ref . jl^] it follows that the 
expansion in A is well defined only if the gauge transformations of the fields are homogeneous 
{i.e. the theory must be gauge invariant order by order in A). This can always be achieved 
via appropriate redefinitions of the collinear fields. In Ref. 3| it has been shown that ho- 
mogeneous soft gauge transformations cannot depend on n-x, hence, in-d derivatives can be 
inserted without spoiling soft gauge invariance. From these considerations it becomes clear 
that matrix elements of the coordinate-space version of our operators will have the form 



{M{pM)\Ux)T'Uym oc 



dw e 



i{upM-x+upM-y) 



(47) 



IV. COMPUTATION OF THE MATRIX ELEMENTS 



In this section we apply the technique introduced above to extract the matrix elements 
of the leading and sub leading SCET operators listed in Sec. |Tl| First we use Eq. fISH)) 
to project the three- and two-particle QCD operators onto SCET. Then we compare the 
resulting expansion with the definition of the light-cone wave functions of pseudo-scalar and 
vector mesons and extract the SCET matrix elements. 



A. SCET decomposition of the QCD operators 

Let us now consider the SCET decomposition of the five independent QCD currents and 
of the relevant three-particle operators, respectively. The matrix elements of the vector, axial 
and tensor currents involve twist-2 and twist-3 light-cone wave functions and it is therefore 
necessary to consider the full expansion in Eq. (j38|) . On the other hand, the matrix elements 
of the scalar, pseudo-scalar and three-particle operators involve only twist-3 and higher wave 
functions and we need to keep only the leading term in Eq. (j2Hl)- 

The SCET decomposition of the various currents are (with z = x — y) 



q{x)q{y) 



(48) 



q{x)-i^q{y) 



(49) 



q{x)-i^'q{y) 



j D2^|jv^(c^)n^ + Vf(cD) 

Ul X^a - ^2 y±a) V^(u;) + (t^i X^a + ^2 y±a) V2"(u;) 



2 



— n^n ■ z 



dt {tx^^ + iy^^) / V'uui V^iu) 



(50) 



10 



Qixh^l q{y) 



z r 



(51) 



(l{x)Of,uq{y) 



-in ■ z{nf,g±^y - n^g^f,^) j dt {tx^a + iy±a) j 'D'^^^ ^ '^'^"(^) ' (^2) 
where g'^ = g^" — }^n^n^ — ^n'^n^, and 

V^UJ = [ dUi duJ2 dUs e5(^i n-x+^,nitx+ty)-^2 n-y) ^ ^53^ 

There are onlythree 3-particle operators whose matrix elements involve twist-3 light-cone 
wave functions l23l. 12511 : 



q{x) gG^^{tx + ty)'jaq{y) = -ina[n^g±^^ - n^g±^y] / "D^a; y V^' (w) , (54) 

(55) 



q{x) gG^uitx + ty) 7^75 q{y) = -i [n^g±^^ - Uyg^^^] / V'^Q ^ -^li^) ^ 



q{x) gGf,u{tx + ty) a^p q{y) = [n^g±^^ - riyg^^^] [riagLpa - n^gxaa] / 'D^uj y T^ito) .(56) 



B. Pseudo-scalar meson matrix elements 

In the following we collect the matrix elements between the vacuum and a pseudo-scalar 
meson state of all the relevant QCD operators [2flll2l|. including the contributions from all 
wave functions of twist 2 and 3: 

(57) 
(58) 
(59) 
(60) 



(61) 
(62) 



{Pip)\q{x)q{y) 


0) = 


= 0, 


{P{p)\q{xh^q{y) 


0) = 


= -ifpf^p / 
^0 


{PipMx)Yqiy) 


0) = 


= 


(P(p)|g(a;)7Vg(l/) 


0) = 


= -^/pp^ 

Jo 


{P{pMx)a^''7,q{y) 


0) = 






{P{pMx)a^''l59G''^{tx + ty)q{y)\0) 

= i[Pf^{Pagu/3 - Ppgau " (/^ ^ '^)]/37r j V^aip^^{d) 
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where fp is the decay constant, /ip = rn?p/{mi + (mj denote the masses of the hght 
valence quarks), /stt is the three-particle decay constant (with dimension GeV^). The 3-body 
integration measure is defined as 



Jv'^a = dai daa dag e'("iP-"+"3P (*"+*^)+"2P J')5(l - ^ 



(63) 



The wave function if{u) is of twist-2 while <fp{u), (fcr{u) and (fsTr{<S) are twist-3. In the 
matrix element of the tensor current one can use partial integration to remove an explicit 
n ■ z factor, which gives 



{PipMx)a>^''qiym 



6 



n ■ p 



,(64) 



Also, the matrix element of the 3-parton operator can be written as 



From the comparison with Eqs. ()48|) - ()56|) we can easily extract the matrix elements of all 
the SCET operators introduced in Sec. |n] 



{p{p)\jA{Qm 
{p{p)\v{Qm 

{p{p)\v{am 

{p{p)\r-^{Qm 

(P(p)|f"^(^)|0) 

{p{p)\r-^{Qm 



-^fpn-p^piu) 
-ifp^p(fp{u) , 



6 



12 



0/3 



—Lp^{u) - <oR 



uu 



u 



u 



—(fcr{u) — 6R 
UU 



G%{u) GZ{u) 



'Py\ 



u 



u 



1 . a/3 V537r(a 



a 



2 ' 



(66) 
(67) 
(68) 

(69) 

(70) 
(71) 



where R = /s^/ (fpfip), {uJi, UJ2) = {u, -u)n ■ p, (cui, uj2, UJ3) = (ai, -02, -a3)n ■ P and 



du 



G%{u) = — I dai I da2 



du 



u — a2 



(72) 
(73) 



All the other matrix elements vanish in the limit in which we neglect twist-4 wave functions. 

Some comments on the manipulations we used to derive the above results are neces- 
sary. The operator V does not appear in Eqs. PH j) - (f^ but is present in the projection of 
q{x)a^^ q{y) ^ hence it can be expressed in terms of the above introduced wave functions. 
The functions G^pJ^u) and G^py{u) have been obtained by inserting Eq. (f7T|) into Eq. (j^ . 
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using the identities 



dtt'' / V^aT{a) 



r 1 

dai / da2 — 

Jo "3 

r 1 

dai / da2 — 
io "3 



u — ai 

"3 
M — 0^2 
"3 



(74) 
(75) 



which are vahd for any function jF(a) and, finally, integrating by parts. The subscripts x 
and y indicate that C'^p^ and Gp^ stem from terms proportional to x_l and y^, respectively. 
The symmetry property of V537r(a) under the exchange ai ^2 implies relations between 



these two functions. For the pion, G-parity implies [21 



V93^(ai, ^2, as) = </'37r(a2, "1, as) 



{*)/ 



U 



(76) 



Finally, the equations of motion in QCD imply relations among the twist-3 wave func- 
tions. In the approach discussed here they follow from the relations among SCET operators 
presented in Sec. II. The first such relation is obtained by noting that the matrix element of 
T^'^ can be directly extracted from the first line in Eq. (j52|) . Comparing this determination 
with Eq. (jnni) we find 



u — U 
uu 



u 



u 



0, 



(77) 



A second equation of motion follows from using the results (|H7j) . (fTHjl in the relation Eq. (fHHjl 
and reads 



(j)p{u) 



1 



Quu 



1 



u u " 



(78) 



Neglecting the terms proportional to the 3-parton wave function 03^, these relations give 
the so-called Wandzura- Wilczek relations for the 2-parton twist-3 wave functions . They 
can be solved exactly and give the solutions (t>p{''j)\ww = 1) (paiu)\ww = Quu. 



Vector meson matrix elements 



Let us collect the matrix elements between the vacuum and a vector meson state of all the 
relevant QCD operators ^ 
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We will work in a reference frame where the momentum of 
the vector meson has a large component along the n direction, and a vanishing transverse 
momentum, thus = ^n^. We neglect a term proportional to the mass squared of the 
meson, because it is irrelevant for the analysis of twist-2 and twist-3 distribution amplitudes. 

The large light-cone momentum component is given by n • p = 2Ey + 0{my / Ey)- The 
polarization vector rj is decomposed as a sum of longitudinal rjw and transverse rjx components 



n ■ rj 



(79) 



Note that these definitions are slightly different from the conventions used in Refs. [23|, |2f 
hence, the expressions for the various matrix elements will look slightly different. Our choice 
of the transverse plane is forced by the structure of the subleading SCET operators. 
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The QCD matrix elements read (we often integrate by parts to remove n • z factors) : 



(80) 



{V{p,r])\q{x)q{y)\0) = —f^ z ■ r)* due'^''P-''+''P-y^ h\^\u) 

2 n-p Jo " 

{V{p,vMxh'q{ym = 0, (82) 



{V{p,vMx)rq{ym = fvmy /'dwe^("^'-^+^^'-^)(r;^*5f)(w) +p^^(0y(M) - («)) 

Jo ^ z ■ p 



fvmv 



-^n^n-pr]*j_ ■ z±J 



_F{u)) , (83) 
{V{p,rjMx)ri'q{ym = -^fvrriye''''^^rj:p,z, jyue'^^^-'^+^P-y^ g'f\u) (84) 

-^n^e'^z^,r^l^g^l\u)\, (85) 







(86) 



(n^n*^ - n'^n^) /if V) - ^m%-r/* {n^'z'i - n'^z"^) / dt;(/i[| V) - <P±{v)) 

(F(p,77)|?(x)^G'^V?(z/)|0) = ^/ymy(ri.p)V(nV*-^'^0 / 2^'5V(c?) , (87) 

{V{p,riMx)gG'^''ri'qiym = \fvmv{n-p)h^^^n-n,rjl, J V'aA{a), (88) 

(y(p,r;)|g(a;)cT"'^^G'^'^g(y)|0) = ^/^m^n-r^* n-pfn"n^^f - n^^n^^^ - n^n'^^ 

8 V 

+nV57) j V^aT{a) , (89) 

where F(u) — Jq dv [(l>\\(v) — g^l\v)] and fv and are the vector meson decay constants, 
defined as 

{V{p, 77)|?7m?|0) = fvmvV*^ , {V{p, v)\qcr^^q\0) = -ifvivlPu - vtPn) (90) 
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In the above formulae, (f)\\{u) and 4>±{u) are of twist-2, while all the other wave functions 
are of twist-3. Since f]±/f]\\ ~ rriv/Ev, we keep terms proportional to my only when they 
involve t]\\. 

The extraction of the SCET matrix elements proceeds in complete analogy with the 
pseudoscalar meson case. Comparing Eqs. (|HT |) -(|S m) with Eqs. (PHjl - lj^Uj) . we obtain: 



{Vip,r^)\Mu)\0) 
{Vip,v)\MCu)\0) 
{V{p,r^)\J^{u)\0) 



-fvmvn-r]*(j)\\{u) 
0, 



(91) 
(92) 
(93) 



(\/(p,r^)|Vr(cD)|0) 
(y(p,r/)|Vr(u;)|0) 

(V^(p,r/)|K(a;)|0) 



-fvmvvT 9±^ (m) 



u — u ^, , 
—F{u) 



uu 



u u 

±F(u) - ^ + ^ 

uu u u 



(94) 
(95) 

(96) 
(97) 



(\/(p,r/)|^^(cD)|0) 
{V{p,r^)\At{Qm 
{V{p,r^)\A'',{Qm 

{v{p,ii)\A',{^m 



\fvmve^lril^g^l^ {u) 
ifvmv e'^r]l^ 9±\u) , 



-fvmv e'l r]l^ 



u-ugt\u) G'Ziu) gP^{u) 
uu 4 u u 

1 9i\u) G^^liu) ^ G^yl{u) 



uu 4 



u 



u 



(98) 
(99) 

(100) 
(101) 



{V{p,^)\S{um 

{v{p,ii)\s{Qm 

{V{p,r^)\VmQ) 

iy{p,r^)\v{Qm 



1 n-p " 

fWv'^h^^{u), 
n-p " 

0, 
0, 



(\/(p,r/)|T^-(cJ)|0) 



{V{p,ii)\T^''{^m 



\fT 2 n-y* uu 
c.Jv''''V - 9± 

2 n-p 



X 



u — u 



uu 



1 J, 2 IT-'V* iiu 

7;fv^v- — 91 

2 n-p 



G^vliu) 



2u 



(102) 

(103) 

(104) 
(105) 

(106) 



2u 



(107) 



X 



uu 



dt;(/i(|^(t;)-0x(t^)) + 



2u 



G%{u) 
2u 
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{Vip,v)\V^iu)\0) 



JvmvVl —jT^ 



an 



a 



2 ' 



t: f i' jJ ^3 



(108) 
(109) 
(110) 



where 



^Vy I 



U) 



dai 



da2 - — ^ (V, A, T) (a) , 



;iii) 

:ii2) 



For the case of the p meson in the hmit of isospin symmetry, G-parity requires that the 
3-parton wave functions satisfy ^(ai,a2,«3) = A{a2,0ii,a^) and {V, T}(ai, 0:2, as) 



— {V, T}(a2, ai, as). This imphes the relations G''yy{u) = —Gy'^{u) and Gyy 

Note that Vi and Ai do not appear in Eqs. (j48p - (j56p but their matrix elements are easily 
obtained via Eq. ()30|). The computation of all the other matrix elements is similar to the 
pseudoscalar meson case. 

The insertion of the matrix elements, Eqs. ()94j ) - ()101|) . into Eq. ()31|) results in the following 
relations between the chiral-even light-cone wave functions: 



u 



4 uu uu 



2g'l\u)- 



u u 



1 .(a), , lu-u , , F{u) _ G^^Uu) - G'^liu) G]::>{u) + G^^>{u) 



■9± [u 



— 9r[u 



2 A uu ' ^ ' uu u u 

The identities Eq. (j32j) imply two other relations for the chiral-odd wave functions 

u — u 



h'P(u) 



uu 



—hu '[u) 

du " uu 



dv{h\f\v) - (f)±{v)) 



U 



G^viiu) 



-G%{u) 
u ^ 



(113) 
(114) 

(115) 
(116) 



In the limit in which we neglect higher Fock state contributions (i.e. where we set to zero 
the right-hand sides of Eqs. ()113|) - (|114D ). we obtain Wandzura-Wilczek like relations between 
twist-2 and twist-3 wave functions [29]. They can be again solved exactly 0,11^ and give 
the twist-3 wave functions in terms of the twist-2 ones (j)\\{u),(f)±{u). For the chiral-even 
structures one finds 



9i\u) 



WW 



WW 



L"'0 

u 



dv— h u 



dv 



V 

dv— — 



V 



(117) 
(118) 
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and for the chiral-odd wave functions 

h^\u) _ = (2m- 1 

h\i\u) = 2 



V 



WW 



u I av h u I av 



(119) 
(120) 



V. REPARAMETERIZATION INVARIANCE CONSTRAINTS 

The soft-colhnear effective theory has an additional symmetry, related to the Lorentz 
invariance of the full theory, which was explicitly broken by defining the effective theory in 
terms of the arbitrary light-cone vectors and fi^. This symmetry manifests itself as an 
invariance under small changes in the light-cone vectors and n^, and is usually called 
reparameterization invariance (RPI). The RPI has been shown to impose rather stringent 
constraints on the form of the effective theory Lagrangian and heavy-light currents ,.10.. 
isl . isl . 3]. For earlier uses of Lorentz invariance in hard scattering processes, see Refs. 




In this Section we show that reparameterization invariance can be used to derive strict 
constraints on the Wilson coefficients of the subleading collinear operators considered in this 
paper. In particular, these constraints fix the coefficients of the three-parton operators of 
the form Xn,u]i^'^iDxW]i^.^Xn,u]2^ terms of coefficients of the leading two-body operators 

of type Xn,ujJ^Xn,ui2- 

The explicit form of the RPI constraints for the subleading collinear operators is in general 
process dependent, and depends on the SCET operators which are allowed in the expansion 
of the physical quantity being described (effective Hamiltonian, current, etc.). Despite this 
diversity, there is one general feature which is common to all these situations: the SCET 
expansion must contain (at least) one additional vector z^, in addition to and n^. This 
vector can be for example f^, the heavy quark velocity in problems involving heavy quark 
decay, or a space-time vector z^, describing the nonlocality of a T-product, as for example in 
hard scattering processes. The presence of the additional vector is required by type III RPI: 
the Wilson coefficients of the SCET operators must depend on RPI invariant arguments, 
and it is impossible to form such combinations from n-V and n ■ alone. The situation is 
different in the presence of an additional vector z, when RPI-III invariant combinations can 
be formed as n ■ zn ■ V and n ■ zn ■ . 

In the following we illustrate the form of the RPI constraints using a few examples. We 
will derive the RPI constraints on the Wilson coefficients in the matching of operators with 
scalar and vector quantum numbers, which will be left completely unspecified. 



A. Scalar operator 

Consider first a scalar QCD chiral-even operator S{z), depending on an arbitrary vector 
Zf^. Keeping terms of 0(1, A), the most general form for this operator in SCET contains the 
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terms (integration over uji is implied on the RHS) 



S{z) = C(^i,^2)(ri-z)Jy(^i,^2) + J]A(^i,^2)(^aVr) (121) 

i=l 

2 



i=l 



This expansion contains the most general structures transforming as a scalar under Lorentz 
transformations, and which are chiral-even and invariant under type-III RPI. As mentioned 
above, the Wilson coefficients C, Di, D, E are functions of the type III invariants n ■ zuji and 
n ■ ZUJ2- In the presence of more than one vector z, the possible number of structures on 
the right-hand side is correspondingly larger. The Wilson coefficients will also depend on 
all possible type III RPI invariant combinations n ■ Zifi ■ 7"^^^ The constraints derived below 
can be straightforwardly generalised to such cases. 

We will show in the following that type-I and -II RPI gives stringent constraints on the 
form of the Wilson coefficients appearing in Eq. p2H) . There are two constraints following 
from type-I RPI 

(I-l): (l + uJi-^ + ^2-^) C(^i, U2) - J2 ^2) - Ie{uj,,uj2, 0) = , (122) 

(1-2): Di{ui,U2) = 0. (123) 
and five other constraints from type-II RPI 

(II-l): C(c^i,c^2)- A(^i,^2) = 0, (124) 

(II-2): (coi^ + u;2^] C{u;,, ^2) - A(^i, ^2) = , (125) 



" duji duj2 

(II-3): bi{ui,U2)=0, (126) 
(II-4): (uj^^-uj2^\ C{uj^,uj2)-D2{ui,U2)=Q, (127) 



\ " duJi duj2 J 

(II-5): C(cJi -CU3,CU2) -C(cui,cu2 + cj3) (128) 

Ar\^E{ijJx,^2,^z) =0. 

In the following we present the derivation of these constraints, focusing on the type-II RPI 
constraints, which are technically more involved. We start by computing the variation of the 
leading 0(A°) term in 1)1211) under a type II reparameterization transformation fiy^^e^. 
This receives contributions from the change in the Wilson coefficient and in the operator 

bii{C{n-zUi)Jy{yji)\ = 6u[C{n-zuJi)]Jv{uJi) + C{n-zuJi)6n[Jv{o^i)] ■ (129) 

The variation of the leading order operator J'v{uJi) in Eq. (|129|) can be computed using the 
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action of the RP transformation on the colhnear fields given in [16|. One finds 



1 



n 



2 n-V 



(130) 



n-V 



r 



Xn,t 



= ^£±-Vl{uJi) - j duj3{e±-V3iuJi+UJ3,U2,UJ3) - £^-V3iuJi,UJ2-UJ3,UJ3)} , 

where the first hne comes from the change in the colhnear quark fields, and the second from 
the variation of W. 

The variation of the Wilson coefficient in Eq. ()129|) can be found by recalling that ui — > 
n ■ and uj2 ^ n ■ V 

2' 



dcoi 



;i3i) 



The two terms can further be written in terms of the operators introduced in Sec. [H] as 

(132) 

duJ3UJ3e±-V3{uJi + UJ3,UJ2,UJ3) 



-uJi[e±-V2 + e_L-V2] 



and 



Xi^i^(^±-'P±)Xu 



Aw^e^-zD,^^^, - x^,^[W^e^-{iD,^yC,U (133) 



Xi 



-UJ2[e±-V2 - e±-V2] - I duJ3UJ3e±-V3{uJi,UJ2 - UJ3,UJ3) . 



Combining everything, one finds the following total result for the variation of the leading 
term in Eq. ()12H1 under a type-II RPI 



(134) 



d 



-UJl[e±-V2+ e±-V2]- I duJ3UJ3e±-V3iuJl+UJ3,U2,UJ3) 



+ 



d 

dU2 



1 



duJ3^3e ^■'V3{yJ\-,^'l -UJ3,^3) 



£:_L-Vi(t^i) - / du)3{e^-V3{uJx\ uj3,uj2,uj-^ - e^-V3{uJx,^'i- ^■i.^'i)) 



+C(cJi 

This has to be cancelled by the 0(A) terms in the variation of the power suppressed terms 
in Eq. ()121|) . The transformation of all colhnear subleading operators Vi,2; ^\,2 and V3 is the 
same, and is given by 

(a;,) = -ln"£^-V.(a;.) . (135) 
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Using this relation in Eq. p2ip . and requiring the cancellation of all independent structures, 
gives the type-II RPI constraints Eqs. (IHlll - dT^ . 

The type-I RPI constraints Eqs. p22|) - ()123|) are derived in a very similar way, so we only 
sketch the relevant steps. Under a type I RPI transformation ^ + A^, the variation 
of the leading order term in Eq. (jl2H) is given by 



OUJi OUJ2 



Muj,) . (136) 



This has to cancel the 0(A) term in the variation of the subleading terms. The action of 
the type I RPI transformation on the individual operators is easily obtained as 

5f (c.,) = 5f ^2^(0;.) = -AiMui) , (137) 
5f ^Vf (cu,) = -te^'A^JAiui) , (138) 
6\'^V^{u,)=0, (139) 

6\'^V^{UJ,UJ2UJ3) = -^A>'^Jv{uj,U2)6M . (140) 

Inserting this into Eq. ()12H) . it is easy to see that the total variation of the Vf terms cannot 
be cancelled by anything else. This gives Di = 0. Requiring the cancellation of the terms 
proportional to {z±-A-^)^v gives Eq. (jl42p . This completes the proof of the RPI constraints. 

We note here the remarkable result that the Wilson coefficients of the 0(A) operators 
in the expansion of a scalar operator are completely fixed by RPI in terms of the Wilson 
coefficient C{uJi) of the leading operator. This is different from the case for the heavy-light 
currents studied in Ref. where several of the subleading operators were not constrained 
from symmetry arguments alone, and their Wilson coefficients have to be determined by 
explicit matching computations. 



B. Vector operator 

We consider next the reparameterization invariance constraints on the Wilson coefficients 
in the matching of a chiral-even operator with the quantum numbers of the vector current 
Vn{z). For this case, the SCET expansion can contain more terms. Keeping again terms up 
to subleading order we can write again the most general type-Ill RP invariant structure as 

2 2 



Vf,{z) = n^C{uJi,uj2)Jv{uJi,(^2) + '^Bi{uJi,uj2)Vi^ + '^Bi{uJi,uj2)Vif, (141) 

j=l i=l 
2 2 

n. ■ 7. ^ n. ■ 7. ^ 



n-z — ' n-z 

i=l 1=1 

n 

+ Ei{uJi,UJ2,UJ^)V-ii, H -E2{uJi,UJ2,UJ'i){z-Vz) ■ 

n-z 



We neglect here terms of the form Zfj_S{z) with S{z) a scalar quantity which is RP invariant 
by itself. The most general form for such a term has been discussed in the previous section. 
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The constraints for this case can be derived in analogy with those for the scalar current. 
We find four constraints from type I RPI 

(I-l): C(^i,^2)- 5^ 072) -^^1(^1, 0^2,0) =0, (142) 

i=l,2 

(1-2): (^c^i^ + c^2^^C^(cui,cu2)- ^ A(cui,o;2)-^^2(a;i,c^2,0) = 0, (143) 

(1-3): B,{u,,uJ2)=0, (144) 
(1-4): D,{u;,,u;2) = 0. (145) 

and four other constraints from type II RPI 

(II-l): CiuJi,uj2)-BiiuJi,co2)-DiiuJi,co2)=0, (146) 
f d d \ 

(II-2): [UJI—+LJ2T^]C{UJI,UJ2)-B2{UJI,UJ2)-D2{UJI,U2)=0, (147) 



" duji duj2 

(II-3): (^a;i^-a;2^^ C(a;i,^2)-52(^i,^2)-£>2(^i,^2) = 0, (148) 

(II-4): C{uJi-uj^,uj2) -C{uJi,uj2 + uj^) (149) 
f d d 



+ \ ^ Ei{uJi,UJ2,UJ3) = 0. 



2 

i=l,2 



It is instructive to compare the RPI constraints Eqs. ()142|) - ()149p with the explicit results 
for the matching of the nonlocal vector operator f^(x)7^g(y) onto SCET operators given in 
Eq. (jSni)- For simplicity we take = —y^ = z^j, in Eq. (|^. which can then be written in a 
form similar to Eq. ()14H) provided that one takes 



B,{iu,,iU2) = Ciiu,), 52(a;,) = 0, 5i,2(^^) = 0, (151) 
Di{uJi,uj2) = 0, D2{uji) = -n-z{ui + U2)C{uJi,LJ2) , (152) 

Di{uJi,uJ2) = , D2{uJi) = -n-z{ui-uj2)C{uJi,uJ2) , (153) 



E^{uj^,uj2,uJs) = 0, E2{uJi) = - rft(2t-l)cu32(n-;z)2e^"-^('^i+^2+-3(2t-i))^ ^^54) 

Jo 

Performing the integration in E2{uJi) one finds 

E2{uJi) = -2 ( 1 - -n ■ ;zcu3 ) et"-^(^i+^2+^3) + 2 ( 1 + -n ■ ^cug ) et"-^(^i+^2-^3) , (155) 



It is easy to see that the constraints Eqs. p42|l - p49|l are indeed satisfied with these 
results for the Wilson coefficients. For the case of the nonlocal vector current, the matching 
Eq. can be worked out using the equations of motion in QCD and is exact to all orders. 
The true power of the RPI constraints Eqs. ()142p - p49|) becomes apparent in those cases 
where the Wilson coefficients have to be obtained in perturbation theory, and where such 
constraints can provide a useful check. 
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-P2 

a) 






FIG. 1: QCD diagrams contributing to photon emission from a qq pair created by the weak vertex. 
The cross indicates an insertion of the weak current 7^ . 



Pi 
if 

N 

-P2 N 



(g).y\/w\/w 



b) 



(g)v/V\/V\/W 



d) 



FIG. 2: SCET diagrams. The cross indicates the insertion of ah the leading and subleading effective 
theory operators. The leading and subleading couplings of collinear quarks and gluons are given, 



for instance, in Ref. 11 1. 



VI. APPLICATION: WEAK ANNIHILATION IN 5 ^ ^7 DECAYS 

In this Section we discuss an application where the subleading operators constructed in 
Sec. IIVI appear in a situation of physical interest. Consider real photon emission from a 
energetic qq pair, which subsequently hadronizes into a meson. This process is relevant for 
the rare radiative decay B V'j, where it contributes to the weak annihilation amplitude. 

In the following we will match the amplitudes contributing to this process onto SCET 
operators, including the three-parton operators constructed in Sec. IHl The matrix elements 
of these operators are then computed with the help of the relations in Sec. IIVI For the case 
of a vector current, the corresponding matrix element is fixed by an exact Ward identity 
31 . 3^ . In addition to serving as illustration for the use of the matrix elements computed 
here, this will provide, at the same time, a strong check for the consistency of our results. 

Let us consider the matrix element of the weak current (I'j'^Plu {Pl = (1 — 7^)/2) between 
the vacuum and a state with a transverse polarised vector meson and a photon: 



(l^(p,77)7(g,e)|(rfVPLu)(0)|0) 



—te 



dS e^«-^ r/) I T [( J7'^Pz.«) (0) , 3!^.. 



(156) 



The conservation of the electromagnetic and the weak currents implies in the usual way the 
following Ward identities for the matrix element of the time-ordered product in Eq. (jl56|) |,3ll 



22 



qo.{V{p,r^)\TnO) 



^{Qd-Qu)fvmvr]'"\ 
^{Qd-Qu)fvmvv"* ■ 



(157) 
(158) 



We checked that the electromagnetic Ward identity, Eq. p57|) . is satisfied at leading and 
subleading order but its check does not require any cancellation between the matrix elements 
of the various subleading SCET operators. We will therefore focus onto the weak current 
Ward identity, Eq. (fTH^ . 

In the following we project the T-product ()156p onto SCET operators, and then use the 
matrix elements computed in the previous section to reproduce Eq. p58|) . For generality we 
keep the Dirac structure of the weak current completely general dTu and define the operator 
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For definiteness, we take the photon momentum along the n direction = E^n^, and the 
meson moving with a large momentum component along the direction. Expanding the 
QCD graphs contributing to the T— product and keeping terms up to subleading order gives 
the following tree level matching 



rpa 
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^(1) 



^ (0) + ^ (1) . 
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n-quJi 



(160) 
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+ iQu 
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n-qLU2{uj2 + uJs 



id) Y 



Xn,w\ 



n-V 



/ A7i,a;2 



1^3 



On the right-hand side, integration over uoi is implied. 

Some details about this matching computation are perhaps in order. The two-body 
subleading operators (the first and third lines of T^^-^) can be obtained from expanding the 
graphs with an external du quark pair (see Fig. 1(a)). The remaining three-body operators 
(the second and fourth lines of T^^-^) are computed by expanding the QCD graphs in Fig. 1(b)- 
(d) and subtracting the insertion of the two-body operators computed in the previous step 
(Fig. 2(b)-(d)). 

The case of the weak current considered in Eq. ()156p is obtained by taking F — » ^^Pl and 
expressing the operators in Eqs. ()160p in terms of the subleading collinear operators intro- 
duced in Sec.|n] We will assume the a Lorentz index to be strictly orthogonal, corresponding 
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to real photon emission. This gives 

^ - ^(0) + ^(1) ' 

— + — + — 

Ol UJ2 J \ 
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(163) 
(164) 



,(165) 



After the insertion of the SCET matrix elements and the contraction with p^ + q^^ = + 
E^n'^, the matrix element of the leading order term vanishes, {p^ + q^){Tl^^) = 0. The 
subleading contribution, on the other hand, gives 
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where the integrations du and doi do2, respectively, are understood. Using the 
QCD equations of motion [2^ 2^ it is easy to verify that the e'^i]l/3 pieces vanish while the 
integral over the rj"* terms reproduces exactly the Ward identity, Eq. ()158j) . In proving this 
cancellation, the following identities are useful 
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'oi(l — O2J 



0. 



Adopting the parameterization of the vector meson wave functions given in Refs. in 113, 
we can get an idea of the numerical impact of the three-particle contributions. First of 
all, note that the equations of motion of QCD allow to express the twist-3 two-particle 
wave functions in terms of the twist-2 two-particle and twist-3 three-particle ones. This 
implies that, in the limit in which we consistently set to zero the functions V and A, the 
Ward identity will be still satisfied. This is a peculiar property of the check that we are 
performing and it is not a general feature of exclusive decay amplitudes. From the inspection 
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FIG. 3: Matrix element (modulo a factor fviT^vVi") of the subleading operator V2 with (solid 
line) and without (dashed line) the inclusion of the three-particle distribution amplitudes. Note 
that in the latter case we consistently set to zero all the twist-3 three-particle wave functions. The 
parameterization of the wave functions as well as the numerical value of the various parameters 
correspond to the p meson and are taken from Ref. 24, 2^ (see also Add. IB]l. 



of Eq. ()166p we see that three-particle hght-cone wave functions can enter the amphtude 
either by modifying the SCET matrix elements of the subleading two-particle operators 
{e.g. V2, A2, ...) or with new genuine contributions {e.g. V^, ^3, ...). In our case 
both contributions are of order 10% and cancel each other. We expect, therefore, effects 
of order 20% in generic exclusive amplitudes that start at subleading order in A. This is 
quite common for decays with vector mesons in the final state because, in this case, the 
usual V — A weak current has no leading projections onto states with a transverse polarised 
meson. For pseudo-scalar mesons, on the other hand, we do not expect large contributions 
because the leading SCET operators have non vanishing projection and are expected to give 
the dominant contribution. 

Finally, note that the importance of these contributions strongly depends on the form of 
the hard scattering kernel. Consider, for instance, the operator V2 • We plot in Fig. El its 
matrix element, Eq. (jUUI) . with and without the inclusion of the three-particle contributions. 
In the former case, the shape of the distribution amplitude is changed and configurations 
closer to the center-point {u ~ 0.5) become more important. 



VII. CONCLUSIONS 

In this paper we classify the leading and subleading SCET operators (in powers of A) 
necessary for the analysis of exclusive processes that involve fast moving light mesons. Some 
aspects of SCET at subleading order still have to be fully understood, but, in order to 
perform our analysis, one only needs the field content of the theory and the corresponding 
gauge transformations. At subleading order, both two-particle (~ ^'^d three-particle 
(~ ^A^O operators are present. 

We show how to express the matrix elements of all the SCET operators between the 
vacuum and a meson state in terms of the standard two- and three-particle meson light- 
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cone distribution amplitudes. Exact operator identities between various SCET operators 
result in corresponding identities between different distribution amplitudes. We checked 
that all these relations are exactly satisfied by the application of the QCD equations of 
motion. In the limit in which the contribution of the three-particle operators is neglected, we 
recover several Wandzura-Wilczek relations between wave functions of geometrical twist-2 
and twist-3. This reflects the fact that operators with a given A scaling have fixed dynamical 
twist {i.e. they have a fixed A/Q suppression); due to the mismatch between dynamical and 
geometric twist, their matrix elements contain contributions from wave functions of different 
geometrical twist. 

The subleading coUinear structures considered here appear as building blocks for SCET 
operators contributing to power suppressed processes. We showed that Lorentz invariance 
(manifested in the SCET as reparameterization invariance) severely constrains the Wilson 
coefficients of such subleading operators, and connects them to the coefficients of the leading 
operators. 

In order to further validate our results and get a feeling for the size of the three-particle 
contributions, we checked a Ward identity for the final state emission of a photon from a 
transverse polarized vector meson. Comparing the exact result with the computation in 
the Wandzura-Wilczek limit (in which the three-particle distribution amplitudes are set to 
zero), we find that corrections of order 0(20%) at the amplitude level are possible. 

In general three-particle light-cone wave functions are expected to be more important 
in processes involving vector mesons in the final state {e.g. B —>■ K*n, ...). In fact, the 
amplitude for the production of a transverse polarised vector meson does not receive twist-2 
contributions (in SCET language, the matrix elements of the leading order operator are 
vanishing). On the contrary, decays that involve only pseudo-scalar mesons {e.g. B tttt, 
...) start at twist-2 and, therefore, three-particle corrections are small. 
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APPENDIX A: NOTATION AND CONVENTIONS 

We use the standard convention for the metric tensor in Minkowski space gi^^ — 
diag(l, — 1, — 1, — 1) and the totally antisymmetric tensor £^^^f"^ is defined with £0123 = 
_^oi23 _ Yoi the Dirac matrices we use the standard conventions, in particular 
^5 — ^^o^i^2^3 define a^u by 

From this we get the relation 0"^,y7^ = —^e^ypacr^", and its inverse — —■^SaPiiv'^^^lb- The 

dual gluon field strength tensor is defined by Gfj,u = \epupcrG'^" ■ 

The hadrons are taken to be moving with large momentum along the light-cone direction 
n. The opposite direction light cone vector n is chosen such that — n^ — {) and n-n — 2. 
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Momenta are decomposed along the light cone as j9 = {p-,p'^,p±) with 



,^ n-p „ n-p „ 
= -yn^ + + Pi- ■ 

We thus define the perpendicular component of the 7 matrices by 

' 2 2 

Similarly we write the perpendicular component of the metric tensor: 

y± y 2 2 
We define the perpendicular components of the total antisymmetric tensor as 



(A2) 



(A3) 



It satisfies the following relations which are useful in practical calculations 
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2 T±7± = -^91 - 2 7 
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APPENDIX B: WAVEFUNCTION PARAMETERIZATIONS 

We collect in this Appendix simple parameterizations of the light-cone wave functions 
used in the numerical evaluations in the text, following [3, H^. For more details about 
their derivation and the com put ation of the noriperturbative parameters appearing in these 
formulas, we refer to HI, ElM H H H H S, Q . 

We start by listing the pseudoscalar wave functions (f>p{u), 4>a{.u), 4>3n{c(i, a^-, as). Working 
at next-to-leading order in the conformal expansion, they can be written as (with 1^ = 2m — 1) 



0p(u) 
0^(m) 

037r(ai) 



1 + 15i?[3^^ - 1] , 

6M(l-M)|l + (5-^^i,o)i?[ye'-| 



360aia2a3 ( 1 + ^t^i,o(7a3 - 3) 



(Bl) 
(B2) 
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For the vector meson case, we give only the wave functions corresponding to the p meson. 
Working again at next-to-leading order in the conformal expansion, the twist-2 two-parton 
wave functions are 



6mm 



1 + a. 



2 -^Wn 



1] 



(B4) 
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The remaining twist-3 two-parton wave functions g^l'"'\u) and h^^'^\u) can be expressed 
with the help of the equation of motion in terms of the twist-3 three-parton wave functions. 
The leading terms in their conformal expansions are 



{V,T}(ai,a2,a3) 



^(ai,a2,a;3) = 360C3aia2tt3[l + ^(7a3 - 3)] • 



(B5) 
(B6) 



We follow here the notations of j2J, |25[. The numerical values of the nonperturbative 



constants are a\ 



0.18 ± 0.10, 



0.2±0.1,C3 



0.032,^^3^ 



3.8, u;J = 7.0, cu^ 



-2.1 



(at the scale = \ GeV^) [2J|. Using the equations of motion, the chiral-even two-particle 
twist-3 wave functions are determined as 
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The corresponding results for the chiral-odd two-particle twist-3 wave functions are 
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